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ABSTRACT 

We report the discovery of an optical counterpart to a quiescent neutron star in the globular cluster 
uj Centauri (NGC 5139). The star was found as part of our wide-field imaging study of uj Cen using 
the Advanced Camera for Surveys (ACS) on Hubble Space Telescope. Its magnitude and color (i?625 
~ 25.2, i?435 — Re25— 1-5) place it more than 1.5 magnitudes to the blue side of the main sequence. 
Through an Ha filter it is ~1.3 magnitudes brighter than cluster stars of comparable i?625 magnitude. 
The blue color and Ha excess suggest the presence of an accretion disk, implying that the neutron star 
is accreting from a binary companion and is thus a quiescent low-mass X-ray binary. If the companion 
is a main-sequence star, then the faint absolute magnitude {Mq25 — 11.6) constrains it to be of very 

low mass {A4^ O.IAMq). The faintness of the disk (Af435 ~13) suggests a very low rate of accretion 
onto the neutron star. We also detect 13 probable white dwarfs and three possible BY Draconis stars 
in the 20" x 20" region analyzed here, suggesting that a large number of white dwarfs and active 
binaries will be observable in the full ACS study. 

Subject headings: globular clusters: individual (NGC 5139) — X-rays: binaries — stars: neutron — 
techniques: photometric — white dwarfs 



1. INTRODUCTION 

Omega Centauri is a prime target for studies of stel- 
lar collisions. It is nearby (D 2± 5 kpc), massive and, 
despite a relatively moderate central density, has one of 
the highest predicted rates of stellar interactions among 
globular clusters, owing to its very large core. Chan- 
nels for compact binary production include exchange en- 
counters, binary-binary collisions, and possibly tidal cap- 
ture (Hut et al. 1992, Fregeau et al. 2003, and references 
therein). Di Stefano and Rappaport (1994) predicted 
that ~100 cataclysmic variables (CVs) formed by tidal 
capture should be in lo Cen at present. Some compact bi- 
naries may also evolve directly from primordial binaries 
in this cluster (Davies 1997). 

We are using the Chandra X-Ray Observatory and 
the Advanced Camera for Surveys (ACS) on the Hub- 
ble Space Telescope (HST) to search for compact binary 
stars in w Cen. The nine fields observed with ACS form a 
mosaic that encompasses more than 100 of the Chandra 
sources identified with ACIS-I in the direction of u Cen 
(Cool, Haggard, & Carlin 2002). The ACS mosaic con- 
tains over a million stars and represents the most com- 
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plete census of stars in u Cen yet obtained. 

Here we report the first detection of an optical coun- 
terpart for a Chandra source using the ACS data. The 
source in question was identified as a possible transient 
neutron star in quiescence by Rutledge et al. (2002) in 
a spectral analysis of 32 of the brightest sources in the 
Chandra data. Its X-ray spectrum was found to be con- 
sistent with thermal emission from a neutron star with 
a hydrogen atmosphere and inconsistent with several 
other possible explanations. The characteristics of the 
optical counterpart that we have identified suggest the 
presence of a disk and binary companion, as in a quies- 
cent low-mass X-ray binary (qLMXB). While more than 
two dozen qLMXBs have been identified in ten different 
globular clusters from their X-ray spectra (Heinke et al. 
2003b), this is only the second optical counterpart found 
for one of these objects during quiescence (the other be- 
ing X5 in 47 Tuc; Edmonds et al. 2002). qLMXBs are 
believed to play an important role in the production of 
millisecond pulsars; studying them in globular clusters 
also holds the promise of new constraints on neutron star 
equations of state (Brown, Bildsten & Rutledge 1998, 
Heinke et al. 2003a). 

We describe the ACS observations and our astromet- 
ric and photometric analyses in § 2. In § 3 we present 
the proposed optical counterpart and discuss it in § 4. 
Results of our search for additional optical counterparts 
to Chandra sources, as well as our study of stellar pop- 
ulations in UJ Cen using the ACS data, will appear in 
subsequent papers. 

2. OBSERVATIONS AND ANALYSIS 

The ACS observations were made on June 27-29, 2002. 
They consist of a 3x3 mosaic of 9 pointings with the 
Wide Field Camera (WFC) covering about 10' xlO', out 
to beyond uj Gen's half-light radius (rh = 288"; Harris 
1996). The X-ray source of interest here lies ~4{5 west of 



2 Haggard et al. 



Table 1. Astrometry 





Chandra R.A. 


Chandra Dccl. 


opt.— X-ray, raw 


opt.— X-ray, corr. 


HST archival 


pixel coords 


Source 


{J2000) 


(J2000) 


Aa, AS 


Aa', A.5' 


image, chip 


X, y 


XA 


13 26 52.141 


-47 29 35.63 


-0'.'08, 0'.'31 


-0'.'03, 0'.'02 


j61p05vsq, 1 


323, 864 


XB 


13 26 53.513 


-47 29 00.37 


-0'.'02, 0'.'28 


0'.'03, -0'.'02 




1005, 1194 


XC 


13 26 48.656 


-47 27 44.88 


O'.'OO, 0'.'33 


0'.'05, 0'.'04 




2682, 280 


qLMXB 


13 26 19.795 


-47 29 10.64 


O'.'ig, 0'.'35 


0'.'24, -0'.'05 


j61p02ffq, 2 


1828, 325 



the cluster center. At each of the pointings we obtained 
3 X 340s exposures with the F625W (i?625) and F435W 
(i3435) fihers and 4 x 440s exposures with the F658N 
(Ha) filter, shifting the camera between exposures to fill 
the chip gap. Having 3 to 4 exposures per filter enables 
us to recognize and eliminate false detections and poor 
measurements caused by cosmic rays. One short expo- 
sure in each of i?625 and B435 was also taken to fill out 
the cluster's horizontal and giant branches. We use these 
filters to look in the X-ray error circle for stars that are 
blue and/or Ha-bright, as potential signatures of accre- 
tion. 

2.1. Astrometry 

Chandra coordinates determined using wavdetect 
fhttpi/Zasc. harvard.edu/ciao) for the quiescent neutron 
star identified by Rutiedge et al. and for three previ- 
ously known X-ray sources in the core, are given in Ta- 
ble 1, columns 2 and 3. To map these positions onto 
the ACS/WFC images, we first applied a distortion cor- 
rection to the individual WFC images using the solution 
obtained for the ^475 filter from a study of 47 Tuc (An- 
derson 2002). This solution should be accurate to ~ 0.15 
WFC pixel (1 WFC pixel = 0'.'05). We then stitched to- 
gether all the individual i?435 frames to make a mosaic 
of the entire field. The images fit together well, with 
no sign of misaligned star images where chips overlap, 
indicating that the distortion correction is working well. 

To determine the R.A. and Dec. associated with stars 
on the mosaic, we used the star lists of Kaluzny et al. 
(1996) and van Leeuwen et al. (2000). More than 15,000 
of the former and 4000 of the latter fall within the field 
of view of our mosaic. There is a small zeropoint offset 
of O'.'S between the two systems, which we take to be in- 
dicative of the uncertainty in the absolute frame; we used 
the Kaluzny system to define the transformation between 
our mosaic system and R.A. and Dec. We estimate that 

the transformation is accurate to ^ 0'.'4 over the entire 
mosaic. 

We cross-checked the optical coordinates using coun- 
terparts previously identified in HST/WFPC2 images for 
two ROSAT sources in the cluster core, XA and XB. We 
recovered stars A and B that Carson, Cool, and Grind- 
lay (2000) identified as CVs, and confirmed that both 
are Ha-bright and blue in the ACS/WFC images (Hag- 
gard et al. 2002). The differences between our optical 
positions for these two stars and the Chandra positions 
(Table 1, column 4) are well within the 0'.'6 uncertainty 
(90% confidence) in Chandra's absolute coordinate sys- 
tem. 

Noting that the X-ray positions for stars A and B are 
offset from their optical positions in nearly the same di- 



rection, we used them to to improve the placement of 
the other X-ray sources on the ACS mosaic. We shifted 
the X-ray positions northwest by (Aa, AS) = (— 0'.'05, 
0'.'30). This shift places star C, identified by Carson et 
al. (2000) as a tentative counterpart of the third ROSAT 
core source, XC, just Ar' = 0'.'06 from the corresponding 
Chandra source position (Table 1, column 5), confirm- 
ing that it is the likely source of the X-rays. Star C is 
Ha-bright but not blue; it may be either a BY Dra-type 
active binary or a cataclysmic variable. 




Fig. 1. — The 20" x 20" region used in the search for the optical 
counterpart to the qLMXB in the -B435 filter. North is up and east 
to the left, approximately. This field is roughly 4^5 west of the 
cluster center. Pre- and post-boresite correction 1" radius X-ray 
error circles are shown (lower and upper circle, respectively). 



The effect of this boresite correction on the position of 
the X-ray error circle for the qLMXB is shown in Fig. 1 . 
The remaining uncertainty in the position of the qLMXB 
is likely to be dominated by the uncertainty in off-axis 
wavdetect positions, ~0'.'5 at the 4!5 off-axis angle of 
the object in the Chandra ACIS-I observation (Feigelson 
et al. 2002). An additional uncertainty of up to '^0'.'4 
associated with the construction of the ACS mosaic and 
the large separation between the qLMXB and stars A 
and B may also be present. We therefore adopted a 1" 
error circle in the search for its optical counterpart. 
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2.2. Photometry 

In view of the variability of the point spread func- 
tion (PSF) in the ACS/WFC (Instrument Science Re- 
port ACS 2003-06), we extracted a 400 x 400 pixel area 
('--^20" X 20") approximately centered on the Chandra 
source position from each of the 12 images (Fig. 1). We 
used DAOPHOT (Stetson 1987) for the analysis, as the 
crowding is significant even outside the cluster core. We 
selected ~10 bright, isolated, unsaturated stars to char- 
acterize the PSF. As it demonstrated little or no vari- 
ability across this small region, we adopted a constant 
model. The best results were obtained by making a sep- 
arate PSF for each image. 

We obtained a preliminary star list using 
DAOPHOT/FIND and then examined and com- 
pared the images by eye in order to remove cosmic rays 
and other spurious detections. We also manually added 
to the star list near neighbors revealed in the course 
of PSF-fitting. Magnitudes and positions for the 1454 
stars identified in this way were then obtained using 
DAOPHOT/ ALLSTAR for each image independently 
Lastly, we matched up stars measured in each of the 
10 long exposures, requiring that the positions match 
to within 1 pixel. We also analyzed the images using 
ALLFRAME (Stetson 1994), which is designed to use 
consistent positions for stars in all images. ALLSTAR 
gave better results, apparently because of the difficulty 
in transforming star positions between non-aligned expo- 
sures with sufficient accuracy, due to the large distortion 
in the WFC. Here we report results obtained using 
ALLSTAR, which we have calibrated using "vegamag" 
zeropoints kindly provided by Sirianni et al. (2003). 
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3. RESULTS 

Color-magnitude diagrams (CMDs) for the 1454 stars 
in the 20" x 20" field are shown in Fig. 2. To reduce 
the effects of any remaining cosmic rays, we have plotted 
median magnitudes. Stars are shown if they appear in all 
3 i?625 and all 3 images and any number of Ha im- 
ages. This ensures that blue stars (e.g., white dwarfs) are 
not overlooked, even if they are too faint to be detected 
in Ha. 

In the B435 — i?625 vs. i?625 diagram (left panel), the 

main sequence can be seen extending ~ 6 magnitudes 
below the turnoff and about a dozen white dwarfs are 

visible at the lower left from Re25 — 23.5 — 25.5 and 
-B435 — i?625 — — 1. In the right panel we plot Ha — i?625 
vs. i?625 • Here the main sequence appears nearly vertical 
and has been shifted so that it is approximately centered 
on Ha — i?625 = 0. In both diagrams, stars within 0''5 
and I'.'O of the boresite-corrected Chandra position are 
indicated with triangles and squares, respectively. 

The most interesting star in the X-ray error circle is 
shown as a filled triangle in Fig. 2. The star is very faint, 
at i?625 = 25.2 and -6435 — -R625 = 1-5, and would not have 
been detected in ^435 or Ha were it not considerably blue 

and Ha-bright. It is ~ 1.5 magnitudes bluer than the 
main sequence and noticeably redder than white dwarfs 

of comparable brightness. In the Ha — R^25 diagram, 
the star lies ~1.3 magnitudes to the left of the main 
sequence, suggesting the presence of a strong emission 
line. Despite being so faint, the star was found in all 
10 images by the automated DAOPHOT/FIND routine, 
making it one of the most reliable detections with i?625 
> 25. No other star in the 1" error circle is significantly 
blue or Ha bright. We identify this star as the probable 
optical counterpart of the X-ray source. 

A finding chart for the proposed optical counterpart 
is shown in Fig. 3. It lies just Ar' = 0"25 from the 
boresite-corrected X-ray position (Table 1, column 5). 
Even without any boresite correction, the X-ray and opti- 
cal positions are offset by 0'.'40 (Table 1, column 4) — well 
within the ~1" uncertainty of HST and Chandra abso- 
lute positions. That the agreement between the optical 
and X-ray positions improves after the boresite correc- 
tion further supports the identification of this star as the 
optical counterpart. 




Fig. 2. - CMDs for stars iu the region shown in Fig. 1. Median 
magnitudes are plotted for stars detected in all 3 i?625 and all 3 
S435 images. Sourees within 0'.'5 and I'.'O of the boresite-eorrected 
X-ray souree position are indicated as triangles and squares, re- 
spectively. The proposed optical counterpart to the qLMXB is 
distinctly blue and Ha-bright (filled triangle). Three BY Dra can- 
didates are also detected (small circles). Approximate magnitude 
limits are indicated with dashed lines. Several probable white 
dwarfs are visible in the lower left of the left-hand panel. 



Fig. 3. — Finding chart for the proposed optical counterpart to 
the qLMXB (see cross hairs) in i?625 (left), B435 (middle), and 
Ha (right) filters. These images are averages of the 3 or 4 available 
exposures in each filter. Error circles of radius 0'.'5 and I'.'O centered 
on the boresite-corrected position are shown. 



We performed several further tests to check on the real- 
ity of this faint object and the reliability of the measured 

magnitudes. First, we examined the star visually in each 
of the ten long-exposure images, and verified that no cos- 



Haggard et al. 



4 

mic rays had affected the photometry. In CMDs made 
from averaged instead of medianed magnitudes, the star 
appears in similar locations. Comparable results were 
also obtained using ALLFRAME, which is somewhat less 
susceptible to contamination of faint star magnitudes by 
bright neighbors. This is reassuring, since the candidate 
has a neighbor just O'.'S (6 pixels) away that is ^^5.5 mag- 
nitudes brighter in i?625 (Fig- 3). We also examined the 
residuals in the vicinity of the star after subtracting it 
out of each of the images and found that the subtractions 
were clean, with no sign of any flux being left behind. 
Thus the object appears to be stellar. 

As a final test, we computed standard deviations for 
magnitudes measured in each filter (Fig. 4). These plots 
give an indication of the accuracy of the photometry as 
a function of magnitude and also show that, in all filters, 

the star is ~ 0.5 magnitudes brighter than the faintest 
stars detected. The accuracy with which it is measured 

(a - 0.2 mag in all filters) implies that unusual -B435 — 
i?625 and Ha — i?625 colors are clearly significant. That 
the star is blue and Ha-bright relative to other stars of 
comparable i?625 magnitude can also be verified visually 
(Fig. 3). 
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Fig. 4. — Plot of sigma vs. median magnitude for the R626) 
-B435 and Ha filters. Detection in all of the "long" exposures in 
a given filter was required for inclusion in the corresponding plot. 
Symbols are as in Fig. 2. Sigma values measured for the proposed 
optical counterpart axe consistent with being due to measurement 
uncertainties alone. 



4. DISCUSSION 

The star we have identified as the optical counterpart 
of the qLMXB is faint, blue and Ha-bright. These char- 
acteristics are typical of semi-detached binaries in which 
a compact object is accreting from a low-mass companion 
(e.g.. Cool et al. 1998). From the observed 1.3 magni- 
tude Ha excess and the widths of the F658N and F675W 
filters, we infer the presence of an emission line with an 
equivalent width of EW(Ha) ~ 220A, after correcting for 



the contribution of the line to the flux through the i?625 
fllter. Such a strong line is reminiscent of lines seen in 
short-orbital-period cataclysmic variables with low ac- 
cretion rates (Patterson 1984, Tylenda 1981). Compa- 
rably strong Ha emission is also seen in the SXT GRO 
J0422-I-32, a black hole candidate with Porb '^5 hrs (Har- 
laftis et al. 1999). Whereas in principle the neutron star 
might have been accreting from the cluster interstellar 
medium (Rutledge et al. 2002), the characteristics of the 
optical counterpart strongly suggest the presence of an 
accretion disk and a low-mass binary companion. 

The object's X-ray to optical flux ratio is high com- 
pared to qLMXBs in the field but is on the order of those 
detected or constrained optically in globular clusters. As- 
suming i\V — R color of ~0.5, we estimate fx/fv ~240 
for this object. Similar estimates for the two qLMXBs 
in 47 Tuc yield fx/fv "^50 for X5 and a lower limit of 
fx/fv ?il80 for X7 (the latter based on a limit on the 
brightness of the optical cotmterpart; Heinke et al. 2003, 
Edmonds et al. 2002). Our limit on the magnitude of the 
optical counterpart of the qLMXB in NGC 6397 implies 
fx/fv >50 (Grindlay et al. 2001). Field qLMXBs Aql X- 
1, Cen X-4 and SAX J1808 have fx/fv -^10 - 20 judging 
from information given by Campana et al. (1998), Cheva- 
lier et al. (1989, 1999), and Homer et al. (2001). 

Further insight can be gained from the intrinsic color 
and absolute magnitude of the optical counterpart. 
Adopting a reddening of E(B-V) = 0.11 (Lub 2002), we 
derive extinction values of A435 = 0.45 and A625 — 0.29 
and thus an intrinsic color of (i?435 — i?625)o = 1-3. This 
color helps to rule out the possibility raised by Rutledge 
et al. that the source could be a narrow-line Seyfert I 
galaxy. For a plausible spectrum, it corresponds to a 
Sloan g' — r' color of ~1.0, which is redder than 150 
such galaxies studied by Williams et al. (2002). The 
high fx/fy ratio further argues against its being an ac- 
tive galaxy, as AGN typically have fx/fv values 1 — 2 
orders of magnitude lower (e.g., Keokemoer et al. 2004), 
as does the excess Ha emission, since a strong line would 
have to be redshifted into the F658N bandpass by chance. 
A measurement of the object's proper motion would help 
confirm its cluster membership, which would also rule out 
the possibility that it is a field qLMXB along the line of 
sight to Lo Cen. We estimate that a baseline of >5 years 
would be required to make such a measurement given the 
faintness of this object and the ~6 mas/yr (~0.1 WFC 
pixels/yr) proper motion of uj Cen. 

As a cluster member, the star is at ^5.0 kpc. Us- 
ing the extinction above, the distance modulus is then 
(m— M)625 = 13.8, yielding an absolute magnitude of 
-^625 = 11.4. Adjusting for the Ha emission line within 
the F625W bandpass gives M625 ~ 11.6 for the contin- 
uum alone. A ±0.5 kpc uncertainty in the distance (cf. 
van Leeuwen et al. 2000 vs. Thompson et al. 2001) intro- 
duces an uncertainty of about ±0.2 mag in these values. 

This is the faintest optical counterpart detected for a 
qLMXB. The only known system that may be fainter 
is the one in NGC 6397, for which a limit of My > 11 
has been derived from its non-detection in HST/WFPC2 
images (Grindlay et al. 2001). The optical counterpart 
to X5 in 47 Tuc and to the field qLMXBs Aql X-1 and 
Cen X-4 (My = 8.2, 8.1, and 7.5 - 8.5, respectively; Ed- 
monds et al. 2002 and references therein) are all consid- 
erably brighter. Such a faint absolute magnitude implies 
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either that the secondary star is of very low mass or that 
it is a compact star. In the case of a main-sequence sec- 
ondary, we can obtain an upper hmit to the mass by 
assuming that it provides all the light in the i?625 band. 
For the metallicity of u> Cen, an absolute magnitude of 
-^625 = 11-6 corresponds to a zero-age main-sequence 
(ZAMS) star with a mass of (Baraffe et al. 

1997). Such a star would contribute about 20% of the 
light in the B435 band, leaving a disk with M435 ~ 13.1. 

We cannot place a lower limit on the mass of the sec- 
ondary star. However, the intrinsic color of the system 
is sufficiently red, despite its location blueward of the 
main seqiience, that it appears likely that at least some, 
if not all, of the light in the i?625 band comes from the 
secondary. To account for its color with a power-law 
spectrum from a disk alone would require a positive ex- 
ponent: f^ ^ A"'^"''. A more typical disk spectrum going 
as would have an intrinsic color of (-B435 — ^625)0 
~0.2 and account for only ~40% of the light in the i?625 
band, even if all the i?435 flux were from the disk. The 
remaining i?625-band light could be accounted for by a 
main-sequence star of mass ^0.12A^q. The difficulty 
of accounting for the object's color without a significant 
contribution from the secondary star argues against the 
possibility that the secondary is compact. 

A main-sequence secondary star mass of ~ 0.12 — 
O.14A^0 imphes a short orbital period for the system. 
In particular, for a ZAMS star in this mass range to fill 

its Roche lobe requires an orbital period of - 1.2 - 1.5 
hr (Warner 1995); longer periods are possible if the star 
is underdense (e.g., Edmonds et al. 2002, Kaluzny & 
Thompson 2003). 

The disk in the system must also be intrinsically very 
faint, with M435 ~13. This is considerably fainter than 
even the disk in SAX J1808.4 — 3658, a transient with a 
2-hr period, whose disk is estimated at Mb = 5.7 — 9.0 
(Homer et al. 2001). Disks this faint are seen in CVs 

with periods below the period gap, i.e., 2 hrs (e.g., 
Sproats et al. 1996), with mass transfer rates of or- 
der of 10~^^ A4q yr~^ driven by gravitational radiation 
(Warner 1995). That an accreting neutron star, with 



a potential well a factor of '^lOOO larger, would have 
such a faint disk suggests an extremely low rate of ac- 
cretion onto the neutron star at present, on the order of 
10~^^ Mq yr~^. This is significantly lower than the rate 
driven by gravitational radiation, if indeed the system 
has a short orbital period. A possible explanation for 
this apparent discrepancy could be that material is be- 
ing transferred to the disk from the secondary, but that 
little, if any, is accreting onto the neutron star from the 
disk at the present time. Such a picture would fit with 
the scenario proposed by Rutledge et al., in which the ob- 
served X-ray emission arises from the neutron star itself, 
and not from present-day accretion. In the Brown et al. 
(1998) model, deep crustal heating of the neutron star's 
core occurs during transient accretion events; a time- 
averaged mass transfer rate of ~8 x 10~^^ Mq yr~^ is 
required to explain the current X-ray luminosity of the 
neutron star in this system. We observe no variability 

in the optical flux (~ 0.2 mag, 1 a — Fig. 4), which is 
also consistent with little or no current accretion onto 
the neutron star. However, the constraints we place are 
limited by the faintness of the star, the small number 
of exposures in each filter and the short (~1.5 hr) time 
span of the observations. 

Finally, we note that several other stars of interest ap- 
pear in the CMD in Fig. 2. Three stars appear mildly 
Ha-bright in the right panel but are on or near the main 
sequence in the left panel (small circles). Visual inspec- 
tions suggest the stars are well measured. These may 
be BY Dra-type active stars as have been seen, e.g., in 
NGC 6397 (Taylor et al. 2001). About a dozen faint blue 

stars with B435 — i?625~ 1 and Rq25 ~ 23.5 are probable 
white dwarfs. Finding such a large number in this small 
region suggests that of order a thousand white dwarfs 
should be detected in the full ACS mosaic. 

We gratefully acknowledge discussions with Scott 
Anderson and thank Marco Sirianni for providing 
ACS/WFC calibration information in advance of pub- 
lication. This work was supported by NASA grant HST- 
GO-9442. 
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